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■ Abstract 

^f) . We study the classical and quantum models of a Friedmann-Robertson- Walker (FRW) cos- 

CN \ mology, coupled to a perfect fluid, in the context of the scalar-metric gravity. Using the Schutz' 

representation for the perfect fluid, we show that, under a particular gauge choice, it may lead 
^ ' to the identification of a time parameter for the corresponding dynamical system. It is shown 

O^' that the evolution of the universe based on the classical cosmology represents a late time power 

^ ' law expansion coming from a big-bang singularity in which the scale factor goes to zero while 

bX), the scalar field blows up. Moreover, this formalism gives rise to a Schrodinger-Wheeler-DeWitt 

(SWD) equation for the quantum-mechanical description of the model under consideration, the 
^S) . eigenfunctions of which can be used to construct the wave function of the universe. We use the 

K*" ' resulting wave function in order to investigate the possibility of the avoidance of classical sin- 

gularities due to quantum effects by means of the many-worlds and ontological interpretation of 
quantum cosmology. 
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1 Introduction 



Classical and semiclassical scalar fields play an essential role in unified theories of interactions and 
' also in all branches of the modern cosmological theories [1]. From a cosmological point of view, 

there is a renewed interest in the scalar-tensor models in which a non-minimal coupling appears 
between the geometry of space-time and a scalar field [2]. This is because that a number of different 
scenarios in cosmology such as spatially flat and accelerated expanding universe at the present time 
[3] , inflation [4] , dark matter and dark energy [5] , and a rich variety of behaviors can be accommodated 
phenomenologically by scalar fields. Traditionally cosmological models of inflation use a single scalar 
field with a canonical kinetic term of the form l/2g^'^ d^(j)dy(j) with some particular self-interaction 
potential V{(j)) like l/2m?(f)'^ or A(/)^, etc. Such a scalar field is often known as minimally coupled 
to the geometry. On the other hand, the scalar field in what is qualified to be called the scalar- 
tensor theory is not simply added to the tensor gravitational field, but comes into play through the 
non-minimal coupling term [6]. 

In this paper we shall study the classical and quantum time evolution of a flat FRW model with 
a perfect fluid matter source and a non-linear self-coupling scalar field minimally coupled to gravity. 
In our model the the scalar field is coupled to the metric through a coupling function F[(j)) in order 
to obtain a scalar-metric theory which as we shall see leads to satisfactory cosmological results. The 
classical cosmology of such a formulation is studied in [7] and it is shown that such scalar-metric 
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coupling gives a satisfactory description of the weak fields and a possible way to remove the missing 
matter problem in non-fiat cosmologies. Also, a generalized version of the model studied in [7] in 
classical and quantum cosmology level is considered in [8]. 

Here, we first consider a gravitational action in the FRW background in which a scalar field 
is coupled to the metric with a generic form of a F{(/)) function. To make the model simple and 
solvable, after some steps, we take a polynomial coupling of the form -F(0) = A^"*. For the matter 
source of gravity, we consider a perfect fluid in Schutz' formalism [9]. The advantage of using this 
formalism in our quantum cosmological model is that, in a natural way, it can offer a time parameter 
in terms of dynamical variables of the perfect fluid [10]. Indeed, as we shall show, after a canonical 
transformation the conjugate momentum associated to one of the variables of the fluid appears linearly 
in the Hamiltonian of the model. Therefore, canonical quantization results in a Schrodinger-Wheeler- 
DeWitt (SWD) equation, in which this matter variable plays the role of time. In terms of this time 
parameter, we shall obtain the dynamical behavior of the cosmic scale factor and the scalar field. We 
show that the evolution of the scale factor represents a late time power law expansion coming from 
a big-bang singularity. Also, the classical behavior of the scalar field shows a blow up in this regime. 
Finally, we deal with the quantization of the model, and by computing the expectation values of the 
scale factor and the scalar field and also their ontological counterparts, we show that the evolution 
of the universe according to the quantum picture is free of classical singularities. 



2 The classical model 

In this section we consider a FRW cosmology in which a scalar field which is coupled to the metric. 
Also, a perfect fluid with which the action of the model is augmented, plays the role of the matter 
part of the model. In the context of the ADM formalism the action reads (in what follows we work 
in units where c = H = IGttG = 1) 

S= f d^Xy/^[R-F{cP)g''''cP,M + '^ I dWhhabK'^'' + I d^x^^p, (1) 

Jm JdM Jm 

where R is the scalar curvature and F{(f)) is an arbitrary function of the scalar field 0. Also, is 
the extrinsic curvature and hab is the induced metric over the three-dimensional spatial hypcrsurface, 
which is the boundary dM of the four-dimensional manifold M. The last term of (1) denotes the 
matter contribution to the total action where p is the pressure of perfect fluid which is linked to its 
energy density by the equation of state 

p = ap. (2) 

In Schutz' formalism the fluid's four-velocity is expressed in terms of five potentials e, (, /3, and S 
as [9] 

C/. = -(e,. + CA- + 05,), (3) 

where ijl is the specific enthalpy, the variable S is the specific entropy while the potentials C and (5 are 
related to torsion and are absent in the FRW models. The variables e and 9 have no clear physical 
interpretation in this formalism. The four-velocity satisfies the condition 

U^U^, = -1. (4) 

We assume that the geometry of spacetime is described by the FRW metric 



ds' = -N\t)dt' + a\t) 



1 — fcr^ 



+ r^(dr + sin^ 



(5) 



where N{t) is the lapse function, a{t) the scale factor and A;=l, and —1 correspond to the closed, 
fiat and open universe respectively. To proceed further, we need an effective Lagrangian for the model 



2 



whose variation with respect to its dynamical variables yields the appropriate equations of motion. 
Therefore, by considering the above action as representing a dynamical system in which the scale 
factor a, scalar field ^ and fluid's potentials play the role of independent dynamical variables, we can 
rewrite the gravitational part of the action (1) as 



Sgrav = j dtCgravia, d, (/),(/)) = j dt I Na^ (^R + ^F{cf))<j)^^ - A 



^ 6 / a k Nd 



where we have introduced the definition of R in terms of a and its derivatives as a constraint. This 

procedure allows us to remove the second order derivatives from action (6). The Lagrange multiplier 
A can be obtained by variation with respect to R, that is, A = Na^. Thus, we obtain the following 
point-like Lagrangian for the gravitational part of the model 

Jl^grav = -J^^d^ + ^^Na + ^F{(j))a^4>^ . (7) 

Also, the matter part of the action (1) becomes Smatt = J d'^xdtNa^p, so the Lagrangian density of 
the fluid is JOmatt = Na^p. Following the thermodynamic description of [9] , the basic thermodynamic 
relations take the form 

p = Po(i + n), M = i + n + -, (8) 

Po 

where po and 11 are the rest-mass density and the specific internal energy of the fluid respectively. 
These quantities together with the temperature of the system r, obey the first law of the thermody- 
namics rdS = dH + pd(l / Po) , which can be rewritten as 



TdS = dll + pd (—) = (1 + n)d [inpo + n)l , (9) 



where we have used the equation of state (2). Therefore, we obtain the following expressions for the 
temperature and the entropy of the fluid 

r = l + n, 5 = lnpo"(l + n)=ln-^Po"- (10) 
Now, we can express po and U as functions of p, and S as 



l + n = ^, Po={^ e--/", (11) 



so that with the help of (8), one can put the equation of state in the form 

On the other hand, normalization of the fluid's four- velocity (3), according to the relation (4) implies 
p = {e + OS) /N. Therefore, using the above constraints and thermodynamical considerations for the 
fluid we find 

Cmatt = (^^^^ + ^S) "^"'^ e-''-. (13) 

Let us now construct the Hamiltonian for our model. The momenta conjugate to each of the above 
variables can be obtained from the definition Pq = In terms of the conjugate momenta the 
Hamiltonian is given by 

H = Hgrav + H^att = dPa + ^P^ + eP, + SPs " (14) 

where C = Cgrav + ^matt- Expression (14) leads to 



H = NH = N 



24 a AF{(j))aP ^ 



(15) 



(6) 
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Now, consider the following canonical transformation which is a generalization of the ones used in 
[11] 



T = -Pse-^Pr^"'^^^ , Pt = Pf +^e^, 



e = e-(a + l)^, 



(16) 



Under this transformation Hamiltonian (15) takes the form 



H = Nn = N 



1 P' 



+ 



24 a 4F{(f))a^ 



Pi - Qka + 



Pt_ 

n3a 



(17) 



We see that the momentum Pt is the only remaining canonical variable associated with the matter 
and appears linearly in the Hamiltonian. The setup for constructing the phase space and writing the 
Lagrangian and Hamiltonian of the model is now complete. 

The classical dynamics is governed by the Hamiltonian equations, that is 

a = {a,H} = -^,^, 

Pa = {Pa,H} = N 

^ = {<P^H} = ^^P^, 

■ Pi F'W 



1 pS 

24# + 4F((/.)o- 



(18) 



Pcl> = {P<j>,H} = Nj^j^, 
f = {T,H} = ^, 
Pt = {Pt,H} = 0. 

We also have the constraint equation "H = 0. Up to this point the cosmological model, in view of the 
concerning issue of time, has been of course under-determined. Before trying to solve these equations 
we must decide on a choice of time in the theory. The under-determinacy problem at the classical level 
may be resolved by using the gauge freedom via fixing the gauge. A glance at the above equations 
shows that choosing the gauge A'' = a^°', we have 



N = a 



3a 



T = t, 



(19) 



which means that variable T may play the role of time in the model. Therefore, the classical equations 
of motion can be rewritten in the gauge N = a^" as follows 



12" 



Pa, 



Pd> „3a-3 



(20) 



A - 1 P'^^l p2„3tt-3 

where we take Pt = Pq = const, from the last equation of (18). The two last equations of the above 
system indicate that the field (p obeys the second order equation of motion 



^ + 2F(^^ + '^'-"^a=°' 



(21) 
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which is equivalent to the conservation law, provided the standard perfect fluid energy-momentum 
tensor is introduced. Equation (21) can easily be integrated to yield 

= Ca6("-i), (22) 

where C is an integration constant. Also, eliminating the momenta from the system (20) results 



3-3q 



(23) 



which is nothing but the constraint equation H = 0. With the help of (22) this equation can be put 
into the form 

Gd^ = Ca^^^-^ - 6A;a^" + Poa^°~\ (24) 



where for the flat case k = 0, its solution reads 

^3Po(l-a)2 



a{t) 



it - 5f - 



C 



3(l-t») 



for a 7^ 1, and 



a{t) = ao exp | \l — 1 \ , 
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(25) 



(26) 



for a = 1, with constants 6 and oq. We may set S 



/8C 



so that a(t = 0) = 0. What remains 



V3Po(l-a) ' 

to be found is an expression for the scalar field 0(t). In the following, we shall consider the case of 
a coupling function in the form F{<p) = A^"*. With this choice for the function -F(0) we are able to 
calculate the time evolution of the scalar field as 



^^(^tanh-i 



3A(l-a) 

exp( tanh"^ 

^ 1 A/3A(l-a) 



/8C 



it -6) 



2 

m + 2 



m / -2, 



(27) 



V3Poil-a) 



/8C 



it -6) 



m 



for a 7^ 1 and 



m = 



^^it-6) 



m+2 



m / -2, 



(28) 



</'oeVT(*-^), 



m 



for a = 1. We see that for a ^ 1, the evolution of the universe based on (25) begins with a big- 
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bang singularity at t = and, for a < 1, follows the power law expansion ait) ~ i3(i-a) a,t late 
time of cosmic evolution while the scalar field has a monotonically increasing behavior coming from 
(j) — >■ — oo, reaches zero and then blows up at a finite time. Also, in the case of a = 1 (stiff matter), 
an exponential solution is obtained for the corresponding cosmology in the chosen time variable. 
In terms of cosmic time r] = J Ndt, solution (26) reads as a(r/) ~ r/^^^ which shows a decelerated 
expansion. This not surprising since for a = 1 there is no violation of the strong energy condition, 
and hence an accelerated expansion cannot be obtained. To understand the relation between the 
big-bang singularity a — > and the blow up singularity (f) — > ±00, we are going to find a classical 
trajectory in configuration space (a, ^), where the time parameter t is eliminated. From (22) and 
(24) one has 



da 



(29) 
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where for the case k = 0, after integration reads 



i^ft^sinh- 



3A(1 - a) 




2 

m+2 



(30) 



Equation (30) describes two branches for which ^ (or a non-zero constant if we add a non-zero 
integration constant to the above relation) if a — >■ oo and (f> — >■ ±oo if a — >■ 0. This means that at the 
late time limit, the scalar field approaches a constant value while its blow up behavior corresponds to 
the big-bang singularity. We shall see in the next section that this classical picture will be modified 
if one takes into account the quantum mechanical considerations in the problem at hand. 



3 The quantum model 

We now focus attention on the study of the quantum cosmology of the model described above. We 
start by writing the Wheeler-DeWitt equation from Hamiltonian (17). A remark about our quantiza- 
tion procedure is that the canonical transformation (16) is applied to the classical Hamiltonian (15), 
resulting in Hamiltonian (17) which we are going to quantize. To make this acceptable, one should 
show that in the quantum theory the two Hamiltonians are connected by some unitary transforma- 
tion, i.e. the transformation (16) is also a quantum canonical transformation. A quantum canonical 
transformation is defined as a change of the phase space variables {q,p) — >■ {q',p') which preserves the 
Dirac bracket [12] 

[q,p] = i = W(.q,p),p{q,p)] ■ (3i) 

Such a transformation is implemented by a function C{q,p) such that 



q\q,p) = CqC-\ p\q,p) = CpC'K 



(32) 



This canonical transformation C, transforms the Hamiltonian as H'{q,p) = CH{q,p)C ^. For our 
case the canonical relations [S, Ps] = [e, Pg] = i yield 



[T,Pt] 



-Pse-^p-^''+^\P,^+^e^ 



Pse-^e^ 



e\Ps 



and 



[e,n] = e-{a + l)PsP-\P, 



(33) 



(34) 



which means that the transformation (16) preserves the Dirac brackets and thus is a quantum canoni- 
cal transformation. Therefore, use of the transformed Hamiltonian (17) for quantization of the model 
is quite reasonable. 

3.1 Schrodinger- Wheeler-DeWitt equation 

Since the lapse function N appears as a Lagrange multiplier in the Hamiltonian (17), we have the 
Hamiltonian constraint 71 = 0. Thus, application of the Dirac quantization procedure demands that 
the quantum states of the universe should be annihilated by the operator version of Ti, that is 



H*(a, (f), T) 



1 P' 



24 a AF{4))a? 



Pt 

6ka H — 5— 



*(a,0,r) = O, 



(35) 



where ^'(a, (f), T) is the wave function of the universe. Choice of the ordering a ^P^ = PaO- ^Pa and 
F{(j))~^P2 = P^F{(j))~^P^ to make the Hamiltonian Hermitian and use of the usual representation 



-idq results in 



-1^ 
da? 



,d_ 
da 



6a~^F(0)- 



-1 



+ Qa~ 



F'{ 



a 



F{(t)f d(f> 



24m 



-3a 



d_ 
df 



^(a,0,T)=O. 



(36) 
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This equation takes the form of a Schrodinger equation id^/dT = H^, in which the Hamiltonian 
operator is Hermitian with the standard inner product 



J(a,cf>) 



We separate the variables in the SWD equation (36) as 

*(a,0,r) = e^^^V(a,<^), 

leading to 



^2 



da? da F{(j)) dcf)'^ F{^Y 8^ 



ip{a,(p) = 0, 



(37) 



(38) 



(39) 



where £^ is a separation constant. The solutions of the above differential equation are separable and 
may be written in the form ip{a, 0) = U{a)V{(p) which yields 



,2 



>(0) = O, 



(40) 



6 fi F'W d I ■■2 



where u is another constant of separation. Upon substituting the relation F[(j)) = into the above 
system, its solutions read in terms of Bessel functions J and Y as 



U{a) = a 



3(l-a) ^ / 3(1-0) ^ 



3(l-c») 



/Me 



3(l-a) 



(41) 



v{4>) = 



d\Jm + l I 

m+2 \ 



$5 V3(m+2) 



d-)Ym,+ l ( 
m+2 \ 



/6A 



m + 2 



^ V3(m+2)' 



where Cj(i = 1,2) and di{i = 1,2) are integration constants. Thus, the eigenfunctions of the SWD 
equation can be written as 



V96^ 3(i-c) 



3(l-a) 



3(1 -a; 



a 2 J 



Z/V6A 



m+l 



+2 V 3(m + 2) 



(42) 



where we have chosen C2 = ^2 = for having well-defined functions in all ranges of variables a 
and (f). We may now write the general solutions to the SWD equations as a superposition of the 
eigenfunctions, that is 



^{x,y,T) 



00 pi 



E=0 Ju=0 



A{E)C{u)-^^E{x,y,T)dEdu, 



(43) 



where A{E) and C{v) are suitable weight functions to construct the wave packets. By using the 
equality [13] 

I e-^^ r^+'Mbr)dr = j^^e--a, (44) 

we can evaluate the integral over E in (43) and simple analytical expression for this integral is found 
if we choose the function A{E) to be a quasi-Gaussian weight factor 



AiE) = 



16 



/ V96E 
3(1 -a)2 (3(1-0) 



exp 



327 



3(1 - a) 



(45) 
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which results in 



Jo 



3(l-c«) 



V96E \ 

-a 2 = 



3(1 - a) 



a 



2\/l-L/2 
3(l-c«) 



exp 



a 



3(l-«) 



47-a(l-a)2r 



(46) 



where 7 is an arbitrary positive constant. Substitution of the above relation into equation (43) leads 
to the following expression for the wave function 



^{a,4>,T) = a 



) 2 exp 



g3(l-a) 



47-z|(l-a)2r 



fc{v)c 
Jo 



2j-i—{l-afT 



2\/l-iy2 . 

J m+1 I 777 — 2 0?Z/. (47) 



^+2 I 3{m + 2) 



To achieve an analytical closed expression for the wave function, we assume that the above superpo- 
sition is taken over such values of u for which one can use the approximation — u"^ ~ 1. Now, by 
using the equality [13] 



(48) 



2m-\-3 n In 

and choosing the weight function C(y) = v "1+2 (1 — v Y' \ we are led to the following expression for 
the wave function ^ 



^'(a,^,T) = Na 

X J 2m+a 



27-i-(l-a)2T 



3(l-a) s(m+2)+l 

(f) 2 exp 



,3(1-") 



47-i|(l-a)2r 



I 3(m + 2) 



6A , m + 2 



(49) 



where s and N are an arbitrary constant and a numerical factor respectively. Now, having the 
above expression for the wave function of the universe, we are going to obtain the predictions for 
the behavior of the dynamical variables in the corresponding cosmological model. To do this, in the 
next subsection, wc shall adopt two approaches to evaluate the classical behavior of the dynamical 
variables in the model which lead to the same results. In the many-worlds interpretation of quantum 
mechanics [14], we can calculate the expectation values of the dynamical variables and, in the realm of 
the ontological interpretation of quantum mechanics [16], one can evaluate the Bohmian trajectories 
for those variables. In figure 1 we have plotted the square of the wave function for typical numerical 
values of the parameters. As this figure shows, at T = 0, the wave function has two dominant peaks 
in the vicinity of some non-zero values of a and 0. This means that the wave function predicts 
the emergence of the universe from a state corresponding to one of its dominant peaks. However, 
the emergence of several peaks in the wave packet may be interpreted as a representation of different 
quantum states that may communicate with each other through tunneling. This means that there are 
different possible universes (states) from which our present universe could have evolved and tunneled 
in the past, from one universe (state) to another. As time progresses, the wave packet begins to 
propagate in the a-direction, its width becoming wider and its peaks moving with a group velocity 
towards the greater values of a while the values of remaining almost constant. The wave packet 



^Onc may have some doubts on this final form for the wave function and the following results due to the assumption 
\/l — v'^ ~ 1. Indeed, sinee the wave function can be a complex function, we may extend the integration domain over 
to < f < 00. In this case with a numerical study of equation (47), we have verified that the general patterns of the 
resulting wave packets follow the behavior shown in figure 1 with a very good approximation. 
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6 



6 




Figure 1: The figures show j\['(a, T)]^, the square of the wave function in four different time parameter T — 
0, 20, 40, 60. The figures are plotted for the numerical values 7 = 1, m = 2, s = 1, A = 1, and we set the equation of 
state parameter a = —1/3. After examining some other values for this parameter, we verify that the general behavior 
of the wave function is repeated. 
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disperses as time passes, the minimum width being attained at T = 0. As in the case of the free 
particle in quantum mechanics, the more locahzed the initial state at T = 0, the more rapidly the 
wave packet disperses. Therefore, the quantum effects make themselves felt only for small enough T 
corresponding to small a, as expected and the wave function predicts that the universe will assume 
states with larger a and an almost constant (j) in its late time evolution. 



3.2 Recovery of the classical solutions 

In general, one of the most important features in quantum cosmology is the recovery of classical 
cosmology from the corresponding quantum model or, in other words, how can the WD wave function 
s predict a classical universe. In this approach, one usually constructs a coherent wavepacket with 
good asymptotic behavior in the minisuperspace, peaking in the vicinity of the classical trajectory. 
On the other hand, in an another approach to show the correlations between classical and quantum 
pattern, following the many- worlds interpretation of quantum mechanics [14], one may calculate the 
time dependence of the expectation value of a dynamical variable q as 



<q>{t) 



< ^'1* > 



Following this approach, we may write the expectation value for the scale factor as 



< a > (T) 



which yields 



< a > (T) 



I— a ' 
' 5— 3a ^ 



27 + 
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(1 - afT'^ 



5127 



3(l-a) 



(50) 



(51) 



(52) 



^ V3-3a; 

It is important to classify the nature of the quantum model as concerns the presence or absence of 
singularities. For the wave function (49), the expectation value (52) of a never vanishes, showing that 
these states are nonsingular. Indeed, the expression (52) for a < 1, represents a bouncing universe 
with no singularity where its late time behavior coincides to the late time behavior of the classical 

2 

solution (25), that is a(t) ~ t^^^. Now we can calculate the dispersion of the wave packet in the 
a-direction which is defined as 

(Aa)2 =< > - < a >2, (53) 

using (49) and (52), we get 



(Aa) 



27 + 



5127 



(1 



a)^r2 



3(l-a) 



(54) 



The result is that the wave packet traveling in the a-direction, spreads as time increases and thus 
its degree of localization is reduced. The width of the wave packet evaluated in (54) agree with the 
discussion in the end of the previous subsection. Indeed, we may interpret the above relation for the 
width of the wave function as the coincidence of the classical trajectories with the quantum ones for 
large values of time. Also, the expectation value for the scalar field reads as 



< > (T) 



with the result 



(55) 



< > (T) 



i-s(m+2) 



/ <^-3(m+2)-l 



i Vex 



const.. 



(56) 
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We see that the expectation value of does not depend on time which is just the behavior predicted 
by the wave function of the SWD equation. This result is comparable with those obtained in [15] 
where a constant expectation value for the dilatonic field in a quantum cosmological model based on 
the string effective action coupled to matter has been obtained. From the classical solutions in the 
previous section, it is clear that this is the classical trajectory obtained from (30) in the limit a — )■ oo. 
Therefore, in view of the behavior of the scale factor and the scalar field, the classical solutions (52) 
and (56) are in complete agreement with the quantum patterns shown in figure 1, and both predict 
a (nonsingular) monotonically increasing evolution for the scale factor and consequently there is an 
almost good correlation between the quantum patterns and classical trajectories. 

The issue of the correlation between classical and quantum schemes may be addressed from another 
point of view. It is known that the results obtained by using the many-world interpretation agree 
with those that can be obtained using the ontological interpretation of quantum mechanics [16]. In 
Bohmian interpretation, the wave function is written as 



(57) 



where and S are some real functions. Substitution of this expression into the SWD equation (36) 
leads to the continuity equation 



2a^ 



,dn dS 



dS 



da da ' da^ '^'^^ 'da F{(f)) d(j) d(j) 
and the modified Hamilton-Jacobi equation 

1 1 /dS\ 



12 dndS 



-:tt- 7^ + 



24a\daJ 4.F{(f))a^ 
in which the quantum potential Q is defined as 



dsy 



-3a 



24o^{i-°) 



0, 



on 



Q = 



1 d'^n 



+ 



1 an 



1 



+ 



F'(0) an 



2Aanda'^ 2Aa'^n da 4a^F{(f))n d(f)'^ Aa^Ficp^Q d<p' 
The real functions fl{a, (p, T) and S{a, (f), T) can be obtained from the wave function (49) as 



0, (58) 



(59) 



(60) 



n = 



s(m+2) + l 
2 



47^ + 2i6(l-«)'r' 



6(l-a) 



47a' 



3(l-a) 



J 



6A 



2m + 3 I „ I , , > 

I 3(m + 2) 



in + 2 
) 2 



z3(l-")(l 



a, 



+ 



3a 



3 1 



arctan 



a 



3(1 



a] 



327 



(61) 



(62) 



8l672 + ^(l-a)4r2 

In this interpretation the classical trajectories, which determine the behavior of the scale factor and 
scalar field are given by 

P^ = ^^ P^ = - (63) 



Using the expressions for and Ps in (20), the equations for the classical trajectories become 



-12a 



9 (1 - afTa 



2-3a 



8167^ + ^(1 



a)4T2' 



and 



2F{ 



.3(i-«);. 



Therefore, after integration we get 

o(T) = ao 



0. 



3(1-0) 



(64) 



(65) 



(66) 
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Figure 2: Up: the figures show the classical scale factor (solid line), the expectation value of the scale factor (small- 
dashed line) and the Bohmian scale factor (large-dashed line). Down: the classical scalar field (solid line) and the 
expectation value or its Bohmian version (dashed line) versus time. The figures are plotted for the numerical values 
7 = 6, Po = 1, C = 2, oo = 1, and we set the equation of state parameters q = —1 (left) and a = —1/3 (right). After 
examining some other values for this parameter, we verify that the general behavior of the curves is repeated. 



and 

(/)(T) = const., (67) 

where ag is a constant of integration. These solutions have the same behavior as the expectation 
values computed in (52) and (56) and like those are free of singularity. Figure 2 shows the behavior 
of the classical scale factor (25), the quantum mechanical expectation value of the scale factor (52) 
and its Bohmian version (66) versus time for some typical numerical values of the parameters. Also, 
the behavior of the classical scalar field (27) and its quantum mechanical expectation value is shown 
in this figure. The origin of the singularity avoidance may be understood by the existence of the 
quantum potential which corrects the classical equations of motion. To get an approximate scheme 
of this issue let us neglect the 0-terms in (60) because of the constant value for from (67). Now, 
inserting the relation (66) in (61), we can find the quantum potential in terms of the scale factor as 

It is obvious from this equation that the quantum effects are important for small values of the scale 
factor and in the limit of the large scale factor can be neglected. Therefore, asymptotically the 
classical behavior is recovered. In this sense we can extract a repulsive force from the quantum 
potential (68) as 

= -^ = 4(3a-l)(3a-2)a-^ (69) 
aa 16 

which may be interpreted as being responsible of the avoidance of singularity. For small values of a 
(near the big-bang singularity), this repulsive force takes a large magnitude and thus prevents the 
scale factor (and then the scalar field) to evolve to the classical singularity (a — )• 0, — )• ±00). 
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4 Conclusions 



In this paper we have studied the classical and quantum dynamics of a scalar-metric cosmological 
model coupled to a perfect fluid in the context of the Schutz' representation. The use of the Schutz' 
formalism for perfect fluid allowed us to introduce the only remaining matter degree of freedom as a 
time parameter in the model. In terms of this time parameter, we have obtained the corresponding 
classical cosmology by evaluating the dynamical behavior of the cosmic scale factor and the scalar 
field. We have seen that the evolution of the universe based on the classical picture represents a late 
time power law expansion coming from a big-bang singularity. We then dealt with the quantization 
of the model in which we saw that the classical singular behavior will be modified. In the quantum 
model, we showed that the SWD equation can be separated and its eigenfunctions can be obtained in 
terms of analytical functions. By an appropriate superposition of the eigenfunctions, we constructed 
the corresponding wave packet. The wave function in this case shows a pattern in which there are 
two possible quantum states from which our present universe could have evolved and tunneled in the 
past from one state to another. The time evolution of this wave packet represents its motion along 
the larger a-direction while the scalar field (j) remains almost constant. As time passes, our results 
indicated that the wave packets disperse and the minimum width being attained at T = 0, which 
means that the quantum effects are important for small enough T, corresponding to small a. The 
avoidance of classical singularities due to quantum effects, and the recovery of the classical dynamics of 
the universe are another important issues of our quantum presentation of the model. These questions 
have been investigated by two different methods. The time evolution of the expectation value of the 
dynamical variables and also their Bohmian counterparts have been evaluated in the spirit of the 
many-worlds and ontological interpretation of quantum cosmology respectively. We verified that a 
bouncing singularity- free universe is obtained in both cases. The use of the ontological interpretation 
has allowed us to understand the origin of the avoidance of singularity by a repulsive force due to the 
existence of the quantum potential. The repulsive nature of this force prevents the universe to reach 
the singularity. 
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